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Abstract The paper proposes an improved quantum associative algorithm 

with distributed query based on model proposed by Ezhov et al. We introduce 
two modifications of the query that optimized data retrieval of correct multi- 
patterns simultaneously for any rate of the number of the recognition pattern 
on the total patterns. Simulation results are given. 
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1 Introduction 

Since the last two decades, there is a growing interest in quantum computing, 
due to the improvement in memory size (such as super dense coding) and the 
speed-up in computing time. Three advantages of quantum computing make 
it possible: (i) the quantum parallelism which is expressed in the principle 
of superposition and provides an advantage in processing huge data sets; (ii) 
the entanglement, the strange quantum phenomenon that links qubits across 
distances and provides the possibility of measuring the state of all qubits in 
a register whose values are interdependent; (iii) the unitary of quantum gates 
which ensure reversibility and therefore overcome energy dissipation. 

Associative memory is an important tool for intelligent control, artificial 
intelligence and pattern recognition. Considering that quantum information is 
information stored as a property of a quantum system e.g., the polarization of 
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a photon, or the spin of an electron, several approaches of quantum associative 
memory for pattern recognition have been proposed. 

Perus et al. worked with quantum-wave Hopfield-like algorithm that has 
been successfully tested in computer simulations of concrete pattern recogni- 
tion applications [1,2, '']. 

The model proposed by Trugenberger makes a link between Artificial Neu- 
ral Network (ANN)-like and logic-gate-based branches of quantum pattern 
recognition. The model is related to the fact that a special case of Hebbian 
memory-storage is equivalent to quantum-implementable NOT XOR gate [ ]. 
However, some critical view had been added indicating that the advantage of 
quantum states as memories for classical information is not clearly demon- 
strated. 

Ventura and Martinez have built a model of quantum associative mem- 
ory where the stored patterns are considered to be basis states of the memory 
quantum state {inclusive method of quantum superposition). Unitary rotation 
operators increase the probability to recall the basis state associated with the 
input pattern [12,9,10]. The retrieval algorithm is based on the well known 
Grover's quantum search algorithm in an unsorted database, which is also 
an amplitude amplification of the desired state [G] . In order to overcome the 
limitation of that model to only solve the completion problem by doing data 
retrieving from noisy data, Ezhov et al. have used an exclusive method of quan- 
tum superposition and Grover's algorithm with distributed query ['>]. However, 
their model still produces non-negligible probability of irrelevant classification. 

Recently, Zhou and Ding have presented a new quantum multi-pattern 
recognition method, based on the improved Grover's algorithm proposed by 
Li and Li [~] , which can recognize multi-pattern simultaneously with the prob- 
ability of 100% [13]. The method introduces a new design scheme of initializing 
quantum state and quantum encoding on the pattern set. However, there is 
an important constraint that the rate of the number of the recognized pattern 
on the total patterns should be over ^. 

This paper suggests, through an oracle operator Im, two modifications 
of the query that can optimize data retrieval of correct multi-patterns simul- 
taneously in the Ezhov's model of quantum associative memory without any 
constraint. In the first modification, Im invert only the probability amplitudes 
of the memory patterns states as in the Ventura's model. In the second modifi- 
cation, Im invert the probability amplitudes of all the states over the average 
amplitude of the query state centered on the m patterns of the learning. 

The main steps of the approach of these models are [10]: 

— Construction of a quantum classification system that approximates the 
function from which set of patterns M {memory states) was drawn; 

— Classification of instances using appropriated Grover's quantum search al- 
gorithm (the time evolution step); 

— Observation of the system. 

The paper is organized as follows: Section 2 briefly present basic ideas of 
Ventura's and Ezhov's models respectively. Section 3 is used to introduce our 



Improved quantum associative memory 



3 



approach and show the comparing results with the above mentioned models. 
In Section 5 we summarize the paper. 

2 Ventura's model and Ezhov's model 

The main purpose of the quantum associative memory build by Ventura and 
Martinez is pattern completion [ ]. That is, it can restore the full pattern 
from partial, but exact, part one. The memory use a storage algorithm and 
Grover's quantum search algorithm for retrieving the patterns. 

Grover's quantum search algorithm can be considered as a rotation of the 
state vectors in two-dimensional Hilbert space generated by the initial and 
target vectors [ ]. The amplitude of the target state increases towards its max- 
imum while the amplitudes of other states decreases after a certain number of 
iterations. 

Each neuron of the memory is a qubit that can be in a state |0), |1) or a 
superposed state a|0)+/3|l), with |ap-|-|/3p — 1. |ap and |/3p are respectively 
the probability of state |0) and state As with a register of n qubits, one 
can compute at the same time all the 2" numbers by using the following 
superposition state 

2"_1 2"-l 

a quantum associative memory can learn or store 2" patterns. In the stan- 
dard Grover's algorithm, Eq. (1) is obtained by applying n times the Walsh- 
Hadamard gate 

W=-^(|0)(fc| + (-l)'=|l)(fc|), fc^{0,l}, (2) 
to the initial state |0). 

2.1 Storage algorithm 

To generate a quantum register (1) in the superposition of only desired states 
from an initial state of n qubits of the network [inclusive method), Ventura 
and Martinez used the storage algorithm that they named algorithm of ini- 
tializing the amplitude distribution of a quantum state. It works in a 
polynomial time and separate the initial state into the already stored patterns 
term and ready to process a new pattern term. The main operator of this 
algorithm is the 2-qubit controlled gate state generation [II, ^] , 

CSP = |0)(0| (8)1+ |1)(1| (g)SP diag(I,SP), 
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for m > p > 1, where m < 2" the number of pattern of length n to be store^, 
each p is associate to a pattern. I denotes the two-dimensional identity matrix. 
The operator (3) separates out the new pattern to be store by assigning to it 
small amplitude so that others operators can't act on it. To do that, we use 
three registers of n, ti — 1 and 2-qubits: 

— \x) = \xi . . . Xn) the register where the m patterns of length n will be 
stored; 

^ I5) = I51 • ■ -ffn-i) a register used like workspace to identify and mark a 
particular state; 

— |c) = \c1C2) a register of two-qubits of control, that is the operator CS^ acts 
when |ci) = 

At the end of the algorithm there is no entanglement between the x-register 
and the two others which are respectively at jO*^""^) = |Oi...O„_i) and 
|0®2) = |00). 

The simplified form of the storage algorithm is; 



Algorithm 1 Simplified form of algorithm 4 (Ref. [ ]) 

1: \ip) = \xi . . .Xn,gi ■ ■ ■ gn-1, C1C2) = |0i . . . On, Oi . . . On-1, 00); {Initialize the register} 

2: for m > p > 1 do 

3: \ip) = FLIP|i/'>; {Generate the state} 

4: = CSP|i/)); {Apply CS^ operator} 

5: iV"} = SAVE|i/)); {Save the state} 

6: end for 

7: |V> = N0T,2lV>; 

8: Observe the system. 



The FLIP operator, 

FLIP = CNDT|'^^^^)CNOtO^^^^.), {1<J< n, z„ ^ (4) 

change the qubits state of the x-register when |ci) — |0) so that they cor- 
respond to the states \P) associated to patterns. The SAVE operator makes 
the state with the smaller amplitude a permanent representation of the pat- 
tern being processed and resets the other to generate a new state for the new 
pattern. At the end of the whole process the system is in the state 

called blank memory in the sense that all possible states have the same prob- 
ability of being recovered upon measurement [n] . The number of steps of the 
storage algorithm is 0{mn) which is optimal because reading each instance 
once cannot be done faster than that. 



1 Generally m < 2". 
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2.2 Retrieving algorithm 

The associative memory proposed by Ventura and Martinez uses for retrieving 
information a modified version of Grover's search algorithm of an unsorted 
database. The original Grover's algorithm has been modified in order to include 
cases where not all possible pattern are represented and where more than one 
target state is to be found. The reader is referred to [ : _'] for more details. It 

should be noted that Grover's algorithm use only 0{\J^) steps to retrieve m 
elements in disordered list of iV = 2" elements, while in classics algorithms 
the best use 0{ — ) steps. 



2.3 Ezhov's model 

As mentioned in the last section, the associative memory proposed by Ventura 
and Martinez can only do completion data. That is bits sequence shown to 
the network should be identical to a part of bits of one of memorized patterns. 
In order to overcome this limitation, Ezhov et al. have introduced a metric 
into the quantum search algorithm in the form of distributed queries [ ]. The 
model is able to retrieve memory states with probability proportional to the 
amplitudes these states have in the query. Their quantum memory can retrieve 
valid stored patterns from a noisy data. 

The model use the exclusion learning approach in which the system is in 
superposition of all the possible states, except the patterns states. If M is the 
set of patterns and m the number of patterns of length n, 

1 

In other words, the exclusion approach for the learning pattern included each 
point not in M with nonzero coefficient while those points in M have zero 
coefficients. 

The distributed query is in the following superposed states 

N-l 

\ReqP} = J2 ^eg^k), (7) 

x=0 

where Req^. obey to binomial distribution 

\\ReqP\\^ = a'^«(P^^)(l - a)«-dH(p.^). (g) 

In equation (8) 

— p marks the state \p) which is referred as the query center; 

— < a < 2 is an arbitrary value that regulates the width of the distribution; 

— the Hamming distance dnip^x) = |p — x| between \x) and \p) is an 
important tool which gives the correlation between input and output; 
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— the amplitudes are such that J2x II^^^SlP — 1- 

The corresponding memory's algorithm is give by Algorithm 2 and the asso- 
ciate Brickman's diagram [ ] by Figure 1. 



Algorithm 2 Quantum associative memory with distributed queries (Rcf. [ ]) 

1: IO1O2 ■ • ■ On) = |0); {Initialize the register} 

2: I'P) = A\0) = ^ ^-x^M 1^^' {Learn the patterns using exclusion approach} 

3: repeat 

4: Apply the operator oracle O to the register; 

5: Apply the operator diffusion T> to the register; 

6: i = i + l; 

7: until i > A 

8: Observe the system. 




(d) 
— -< 




(c) 





Fig. 1 The Brickman's diagram [ 1] of algorithm 2. (a) Qubits are initialize to the state |0). 
(b) The operator A performs the learning of the set M using the exclusion approach and 
gives to all the possibles states the same probability amplitude a^. (c) The operator oracle O 
inverts probability amplitude of all the states over the average amplitude of the query state 
by changing ax to Ox — 2Reqx (^^^^o (^^1^)* '^^^ ■ ('^) "^^^ operator diffusion "D inverts the 
probability amplitude of the states of \<I') over their average amplitude and for the others 
over the value by changing the amplitudes of probability Ox to 2mx mja^^ — Ox- 

(e) After A applications of steps (c) and (d) the system is observed. 
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In the Algorithm 2 or in the associate Brickman's diagram of Figure 1, 

— O is the operator oracle which invert the phase of the query state \ReqP), 

= I-{l-e''')\ReqP){ReqP\, (9) 
0:a,^a,- 2Reql {Reqlfa^j , (10) 

where is the probability amplitude of the state \x). 

— V is the operator diffusion which invert the probability amplitude of the 
states of \^) over their average amplitude and for the others over the value 
0. 

V = {l-e''')\^){^\-I, (11) 

/N-l \ 

I> : H> 2mx ( ^ TO*ax J - Ux- (12) 

\x=0 / 

where is the probability amplitude of a state oi \^). 

— yl is the number of iterations that whilst the maximal value of amplitudes, 
which must be as far as possible nearest to an integer, 

1 27r 
A = T{-+a),T=—,aGN, (13) 

4 LJ 

with the Grover's frequency 



2arcsinE,E= , ^ V ReqP. 



(14) 



x=0,x^M 



Example 1 In order to help clarify, consider a 3-qubits memory where the 
patterns for the learning are |010) = |2) and |100) = |4) and the distributed 
query centered on |011) = |3). For a = \, 

\Req^) = ^|0) + + g|2) + ^|3) + ^|4) + ^|5) + ^|6) + g|7), 

(15) 

„ 1 ^ •? 6 + 6-\/3 , , 

oj = 2arcsin ^^^^ = o.637r ^T = 3.17 =^ for a = 1, yl = 4. (17) 
8^6 ' V / 

The steps 4 and 5 of the Algorithm 2 will be repeated 4 times. 
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The operator oracle (10) is 

= 1-:^ (V3|0) + 3|1) + 3|2) + 3^|3) + |4) + V3\5) + V3\6) + 3|7) 
(\/3(0| + 3(1| + 3(2| + 3\/3(3| + (4| + V3{5\ + V3{6\ + 3(7|) . 

(18) 

As after the learning process the quantum memory is in the state 

|f> = ^(|0) + |l> + |3) + |5) + |6) + |7)), (19) 
the operator diffusion (12) is 

-|6) + |7))((0| + (1| + (3| + (5| + (6| + (7|) 

(20) 
(21) 



■D=l (|0) + |1) + |3) + |5) + |6) + |7)) ((0| + (1| + (3| + (5| + (6| + (7|) - I. 

At the end of the 4 iterations the register is in the state 

= -0.257|0) + 0.031|1) + 0.683|2) + 0.531|3) + 0.228|4) 
-0.257|5) - 0.257|6) + 0.031|7). 

The probability to retrieve the memory states, |2) and |4), is 0.683^ + 0.228^ = 
51.85%. As expected, the memory state closest in Hamming distance to the 
query center state, |2), presented the best probability (46.65%). 




Fig. 3 Probability of correct 
recognition for a set of two 

Fig. 2 Probability evolution with the number of example patterns in Ezhov's 

iterations in Ezhov's model for 3-qubits memory model, 
patterns. The solid line represents the probability 
Pc of a correct recognition and the dotted line the 
probability of an incorrect recognition. 
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Figure 2 shows the probabihty of observing the correct recognition upon sys- 
tem measurement versus Grover's search iterations. The sohd hne represents 
the probabiUty of a correct recognition and the dotted Hne the probabihty 
Pyj of an incorrect recognition. Note that the periodic nature of the algorithm 
clearly appears, and it can be seen that the probability of success Pc is maxi- 
mized after four iterations. At ^ = 4, the ratio 

^ = 1.08, (22) 

that could be considered as the recognition efficiency of the memory patterns, 
shows that the confidence that the recognition is correct fair in the Ezhov's 
model. This clearly appears in the Figure 3 that gives the graphic representa- 
tion of the probabilities of correct and bad recognition of this example. 
It should be pointed out that if the approximation 

N~l 

B::.-=Y. ^ » ^' (23) 

is wrongly use in Example 1, the number of iterations increase to Tl = 9 and 
the probability to retrieve the memory states |2) and |4) is reduce to 36.00% 
and therefore the recognition efficiency is 

^ = 0.56. (24) 

3 Improved quantum associative algorithm with distributed query 

In order to improve the quantum associative memory with distributed query 
such that it optimize the probability of retrieving the learned patterns, even 
for the biggest Hamming distance from the query center, we proposed the 
Algorithm 3 illustrate by the Brickman's diagram of Figure 4, with an operator 
Im- 



Algorithm 3 Improve quantum associative memory with distributed query 

1: IO1O2 . • . On) = |0); {Initialize the register} 

2: = A\0) = 'Ylx^M 1^)' {Learn the patterns using exclusion approach} 

3: Apply the operator oracle O to the register; 

4: Apply the operator diffusion V to the register; 

5: Apply operator Xm to the register; 

6: Apply the operator diffusion V to the register; 

7: repeat 

8: Apply the operator oracle O to the register; 

9: Apply the operator diffusion T) to the register; 
10: i = j + l; 
11: until i> A-2 
12: Observe the system. 



Two cases we be will considered for the operator Im'- 



10 



J.-P. TCHAPET NJAFA et al. 



MESURE 



(i) 



Vl-^) (h) 



(f) 



(e) 



0\^) (g) 
\ 

Y 

/ 



(d) 



(c) 



^|0) ^ I*) 



|o> 



(b) 



(a) 



Fig. 4 Brickman's diagram [ ] of Algorithm 3. The steps here are the same as in the diagram 
of Figure 1, with however two modifications: (e) The new operator Xm acts as the operator 
oracle and (f) the operator diffusion V is apply again before the A — 2 iterations. 



CI; Im invert only the phase of the memory patterns states as in the Ventura's 
model, 

lM^l-{l-en\'p){^\, E (25) 

J-a, if |.T)eM 
I Qx II not. 

Va; G M, the grover operator act as 

VIm\^) = (2|tf')(if|- 1 + 21^X^1)1^) 

= 2\^){W\^) -\^)+2\ip){ip\^) = 1^). 

C2: Im is formally identical to the operator oracle O of Eq. (10), 



/N-l \ 

Im -.a^^ax- 2REQx (^^^ 



ith 



(27) 



(28) 



REQxW ^ ^5]a^(''^)(l-a,)"-'^-(''-), (29) 



where we consider that the distribution have k centers and < at < | is 
an arbitrary value that regulates the width distribution around the center 
b. But in this paper, we will consider that 

\\REQJ\^ - - V a"'"(^^"''(l - a')""''"'''''^^. (30) 

beM 
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where m is the number of patterns for the learning, b is an item of the set 
M of patterns, and we choose the case where a' = ai, is the same for aU 
the patterns. 

It is noteworthy that using a straight forward approach to classification 
and employing Grover's search, Ventura have found that the exclusion method 
exhibit the lowest overall probability of irrelevant classification compared to 
inclusion method [10]. This explains why the exclusion method is use in the 
Algorithm 3. 



4 Simulations and results 

Consider the data of Example 1, the evolution of probabilities with the number 
of iterations are plot in Figures 5 and 7 for Cl-algorithm and Cl-algorithm 
respectively. 




Fig. 6 Probability of cor- 
rect recognition for a set of 

Fig. 5 Probability evolution with the number of example patterns for Cl- 

iterations for Cl-algorithm for 7-qubits memory algorithm 
patterns. The solid line represents the probability 
Pc of a correct recognition and the dotted line the 
probability of an incorrect recognition. 



For Cl-algorithm, 

Im=I- (|2)(2| + |4)(4|), A = 25{T= 11). (31) 

Therefore, the steps 8 and 9 of Algorithm 3 are repeated 23 times. At the end 
of the algorithm, the register is in the state 

IV'^^) = -0.137|0) + 0.0231|1) - 0.876|2) + 0.301|3) - 0.292|4) - 0.137|5) 
-0.137|6) -H0.0231|7). 

(32) 
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result result 



Fig. 7 Probability evolution with the number of Fig. 8 Probability of correct 

iterations for C2-algorithm for 7-qubits memory recognition for a set of two example 

patterns. The solid line represents the probability patterns for C2-algorithm. 
Pc of a correct recognition and the dotted line the 
probability of an incorrect recognition. 



The probability to retrieve the raemory states, |2) and |4), is 0.86266^ + 
0.287552 = 85.23%. The ratio 

A . '-^^ = 5.77, (33) 
0.1477 ' ^ ' 

shows how higher recognition efficiency of Cl-algorithm compared to Ezhov's 
model. 

For C2-algorithm, by choosing a' = 0.1, one find 

|i?i;Q) = 0.285|0) + 0.095|1) + 0.607|2) + 0.202|3) + 0.607|4) 

+0.202|5) +0.285|6) +0.095|7), ^ ' 



Im^I- (0.285|0) + 0.095|1) + 0.607|2) + 0.202|3) + 0.607|4) + 0.202|5) 
+ 0.285|6) + 0.095|7))(0.285(0| + 0.095(1| + 0.607(2| + 0.202(3| 
+ 0.607(4| + 0.202(5| + 0.285(6| + 0.095(7|). 

(35) 

and A — A. Therefore, the steps 8 and 9 of Algorithm 3 are repeated 2 times. 
At the end of the algorithm, the register is in the state 

IV'*) -0.107|0) - 0.024|1) + 0.772|2) + 0.358|3) + 0.477|4) + 0.152|5) 

-0.107|6) - 0.024|7). ^ ' 

The probability of correct recognition of memory states |2) and |4) is 0.772^ + 
0.477^ = 82.69%, which is fairly lower than that of Cl-algorithm. 
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The Tabic 1 which gives the summary of the relevant parameters of Ezhov's, 
CI and C2 methods shows that despite the Cl-algorithm is slower than C2- 
algorithm it leads to a better recognition efBciency of memory patterns. There 
is a significant gap between the recognition efBciency of CI and C2 algorithms 
and that of Ezhov's model. 



Table 1 Relevant parameters of Ezhov's, CI and C2 methods for 3-qubits memory data 
of Example 1. 



Method 


A 




Ezhov's 


4 


1.08 


CI 


25 


5.77 


C2 


4 


4.67 



For a better comparison of CI and C2 algorithms, consider a 7-qubits 
memory where the patterns for the learning are states |23), |59), |61), and 
|110) and the distribution query centered on |60). Figures 9 and 10 show the 
corresponding probabilities of observing the correct recognition upon system 
measurement versus the number of iterations of Grover's search for a = 0.15 
and a = 0.40 respectively. 

Table 2 Relevant parameters of Ezhov's, CI and C2 methods for 7-qubits memory data. 



Method 


a 


a' 


A 


Pc 


Ezhov's 


0.15 




32 


< 10% 


CI 


0.15 




20 


57.04% 


C2 


0.15 


0.10 


10 


51.45% 


C2 


0.15 


0.40 


13 


18.35% 


Ezhov's 


0.40 




21 


<40% 


CI 


0.40 




14 


93.22% 


C2 


0.40 


0.10 


20 


48.37% 


C2 


0.40 


0.40 


12 


54.70% 



It can be seen, as summarize in the Table 2, that 

— the smaller the arbitrary value that regulates the width of the distribution 
a, for the three methods 

— the larger number of iterations A; 

— the lower the recognition efficiency; 

— for a giving value of a in the C2-algorithm, the closer a' 

— the lower number of iterations A; 

— the larger the recognition efficiency; 

— for any value value of a € [0, 1/2], 

— the best recognition efficiency of memory patterns is given by Cl- 
algorithm and the poorest by the standard Ezhov's model; 

— C2-algorithm seem to be faster than Cl-algorithm. 
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CI (a=0.15) 



1^ 



20 25 30 

Iteration 



C2(a=0.15, a'=0.4) 



20 25 30 35 
Iteration 



Fig. 9 Probability evolution with the number of iterations for 7-qubits memory patterns. 
The arbitrary value that regulates the width of the distribution is a = 0.15. 



Ezhov (a=a.4) 




10 12 14 16 18 20 22 
Iteration 

C2(a=0.4. a =0.1) 



10 12 14 16 18 20 22 
Iteration 




Cl(a=0.4) 



10 12 14 16 IS 20 22 
Iteration 

C2(a=0.4, a'=0.4) 

1^- 



10 12 14 16 IS 20 22 
Iteration 



Fig. 10 Probability evolution with the number of iterations for 7-qubits memory patterns. 
The arbitrary value that regulates the width of the distribution is a = 0.40. 



5 Conclusion 

In this paper we have proposed a multi-pattern recognition method with a 
good rate of success based on an improved quantum associative algorithm 
with distributed query for any rate of the number of the recognition pattern 
on the total patterns. We have introduced an operator Im which acts as the 
oracle operator by considering two cases: the case CI where Im invert only the 
phase of the memory patterns states as in the Ventura's model; the case C2 



Improved quantum associative memory 



15 



where Im invert the probabihty amphtudes of all the states over the average 
amplitude of the query state centered on the m patterns of the learning. These 
improvements appeared as factors that increased constructive interferences so 
that the number of iterations is considerably reduced. Simulation results favour 
the Cl-algorithm despite his high number of iterations compare to the C2- 
algorithm. Work is in progress to investigate other possibilities of improvement. 
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A Learning algorithms 

The algorithm 4 use the following gates: 

— The hermitian 1-qubit NOT gate 

NOT=|0)(1| + |1)(0| = (^JJ). (37) 

— Two forms of the hermitian 2-qubit controlled gate CNOT 

CNOT° = diag(NOT, I) and CNOT^ = diag(I, NOT). (38) 

— Four forms of of the hermitian 3-qubit Fredkin's gate F (a controUed-controUed NOT 
gate) 

f''° = diag(HOT,I,I,I), F°l = diag(n, HOT, I, I), (39) 
f1° = diag(I, I, NOT, I), f" = diag(I, I, I, NOT). (40) 

— And Zpj the values of different qubits of patterns, where we take = O®"'. 
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Algorithm 4 Learning the set M by initializing the amplitude distribution 
of a quantum state (Ref. [ ]) 

1: \'4>) = |xi . . .Xn,gi . . .g„_i,ciC2) = |0i . . .On,Oi . . .0„_i,00); {Initialization} 

2: for m > p > 1 do 

3: for 1 < J < n do 

4: if Zpj ^ 2(p+i)j then 

5: |V> = CNOtO^^^^.)|^>; 

6: end if 

7: end for 

8: |^> = CNOt0^^^^j|V>>; 

9: = ¥P\ii); 

10: |^>=F™JV>; 

11: for 3 < A: < n do 

12: \^)=<tL-2,u-i\^y^ 

13: end for 

14: |V> =CNOTl^_^_^^^^|^>; 

15: for n > fc > 3 do 

16: H)=<tL-20.-^H)^ 

17: end for 

18: |^>=F:f///j^>; 

19: end for 

20: |V> =N0Tc2|V>; 

21: Observe the system. 



